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Within the Wilson integral formalism, we discuss the structure of nonperturbative corrections to 
the quark form factor at large momentum transfer analyzing the infrared renormalon and instanton 
effects. We show that the nonperturbative effects determine the initial value for the perturbative 
evolution of the quark form factor and attribute their general structure to the renormalon ambiguities 
of the perturbative series. It is demonstrated that the instanton contributions result in the finite 
renormalization of the next-to-leading perturbative result and numerically are characterized by a 
small factor reflecting the diluteness of the QCD vacuum within the instanton liquid model. 

2 I PACS numbers: 12.38.Lg, ll.lO.Gh 

o ■ 
(N : 

c : 

^ , I. INTRODUCTION 



(N 
(N 



Oh! 



X 



The various aspects of the instanton induced effects in high energy QCD processes had been addressed at the very 



' beginning of the instanton era (see, e.g., Ref. [1]), and this study has been continued in later decade [2]. Recently, 



the interest in them has been revived [3-7] , and the hope of the direct detection of the instanton induced effects has 
^ appeared [8] . One of the important questions in the description of hadronic exclusive processes is the behavior of the 
— , form factors in various energy domains. The present paper is devoted to the analysis of the corrections to the elastic 
' quark form factor at large momentum transfer induced by the infrared (IR) renormalon [9] and instanton effects, 
treating the latter within the framework of the instanton liquid model of the QCD vacuum [10]. 

From the theoretical point of view, the form factor analysis requires a perturbative resummation procedure beyond 
the standard renormalization group techniques, since it exhibits a double-logarithmic behavior. In addition to this, 
^-^ , the resummation techniques developed for this particular case can be applied to the study of many other processes 
which possesses the logarithmic enhancements near the kinematic boundaries. On the other hand, in addition to this 
obvious theoretical use, the computation of the quark form factor has important phenomenological applications. A 
similar resummation approach is also used in study of the near-forward quark-quark scattering and evaluation of the 
Oh! soft Pomeron properties where the nonleading logarithmic contributions are quite important [11]. The application 
^ of this formalism in heavy quark effective theory can be also useful [12]. The quark form factor itself enters into 
the various cross sections of high-energy processes [13]. In particular, this quantity finds the most straightforward 
phenomenological application in the total cross section of the Drell-Yan process in the deep inelastic scattering (DIS) 
scheme, which is proportional to the ratio of the timelike and spacelike form factors [14-16]. In this case, the 
exponentiated quark form factor is expressed in terms of the evolution equation, and can be evaluated, in principle, 
■ to any order in perturbation theory. The analysis shows that this is the high-energy asymptotic behavior that is 
important, and all the logarithmic contributions must be taken into account equally, while the power corrections 
could be neglected [16]. The investigation of the electomagnetic quark form factors in moderate and low-energy 
domains can shed light on the problem of scaling violation in DIS and the structure of constituent quarks [17]. 

The color singlet quark form factor is determined via the elastic scattering amplitude of a quark in electromagnetic 
field 

-^M = Fq [bl - P2)^] U(pi)7^w(p2) , (1) 

where m(pi), v{p2) are the spinors of outgoing and incoming quarks. The kinematics of the process is described in 
terms of the scattering angle x- 

coshx=^ = l + ^ , Q'^~iP2^Pir>0 , pl=pl = m'. (2) 
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It is known that the leading large-Q^ asymptotics of the quark form factor is given by the exponentiation of the 

one-loop term [18] 

i^i^nQ^)=exp(-^CWn^^)+o(<ln^"-^^) , (3) 

where A is an IR cutoff parameter. In general, for a correct consideration of the non-leading asymptotic contributions 
one has to resum all perturbative (such as O (a" In^""^ Q^) , O (a" ln^"~^ Q^) , etc) as well as nonperturbative terms. 
An effective framework for rcsummation of perturbative and nonperturbative contributions is provided by the Wilson 
integral approach [19]. Within this framework, the resummation of all logarithmic terms coming from the soft gluon 
subproccsscs allows us to express the quark form factor (3), in terms of the vacuum average of the gauge invariant 
path ordered Wilson integral [20] 

W{C^) = i-Tr(OlPexp (^zg dx^A^{x)j [0) . (4) 

In Eq. (4) the integration path corresponding to considering process goes along the closed contour C^: the angle 
(cusp) with infinite sides. The gauge field 

4(x) = r°A^(x) , T'' = y , (5) 

belongs to the Lie algebra of the gauge group SU{Nc), while the Wilson loop operator Pol's JdxA{x) jjgg funda- 
mental representation. 

In our recent paper [7] , we applied the Wilson integral formalism to evaluation of the perturbative and nonpertur- 
bative contributions to the color singlet quark form factor at the low normalization point fj. of order of the inverse 
instanton size within the instanton liquid model. In the present work, considering the renormalization group (RG) 
evolution equation we extend the analysis to the limit of large momentum transfers focusing on the asymptotic be- 
havior. We show that the nonperturbative effects determine the initial value for the perturbative evolution, find their 
general structure by analyzing the renormalon ambiguities of the perturbative series, and establish the correspondence 
between them and the instanton induced contribution. 

The paper is organized as follows. In Sec. II we reproduce the known results of the perturbative one-loop calculation, 
and derive the evolution equations taking into account the nonperturbative contribution as the initial value for 
perturbative evolution. In Sec. Ill, we study the consequences of the IR renormalon ambiguities of the perturbative 
series and show how the latter prescribes the form of the nonperturbative corrections to the asymptotic behavior of 
the form factor at large Q^. In Sec. IV, these nonperturbative effects are estimated in the weak-field approximation 
within the instanton model of the QCD vacuum. Finally, the large behavior of the form factor is analyzed taking 
into account the leading perturbative, IR renormalon, and instanton induced contributions. The latter are found to 
be determined by small factor expressed via the parameters of the instanton liquid model. 



II. ANALYSIS OF THE PERTURBATIVE CONTRIBUTIONS TO THE WILSON INTEGRAL 



The Wilson integral (4) can be presented as a series 

w{c^) = i + ^{o\Y,{zgr [ I ... [ dx;^ dx;:\...dxl. 



^(x",x"-\...,^i) Tr 4„(x")4„_,(;r"-i)...4,(xi) JO), 



(6) 



where the function 9{x) orders the color matrices along the integration contour. In the present work, we restrict 
ourselves with the study of the leading order (one loop — for the perturbative gauge field and weak-field limit for the 
instanton) terms of the expansion (6) which are given by the expression 



j dxn / dy^ D^^{x - y) 



(7) 
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where the gauge field propagator D^v{z) in n-dimensional space-time (n = 4 — 2e) can be presented in the form 

D^^z) = Ai(£,^2^mVa') - d^d.A2{s,z\iJ,yX^) . (8) 

The exponentiation theorem for non-Abelian path-ordered Wilson integrals [21, 22] allows us to express (to one- loop 
accuracy) the Wilson integral (4) as the exponentiated one- loop term of the series (6): 

W^(CJ=exp[wW(C^) + 0(a^)] . (9) 

In general, the expression (7) contains ultraviolet (UV) and IR divergences, that can be multiplicatively renormalized 
in a consistent way [23]. In contrast to the previous paper [7], we use the dimensional regularization in order to 
work with UV singularities, and define the "gluon mass" as the IR regulator and the parameter /z^ as the UV 
normalization point. The dimensionally regularized formula for the leading order (LO) term (7) can be written as [7] 

W^«(Cx;e,/xVA2,a,) = 87ra,CFMx)(l-£)Ai(e,0,MVA') , (10) 
where h{x) is the universal cusp factor, 

Kx) = Xcothx - 1 , (11) 

and, in case of the perturbative field, 

i,(,,„,,VAn = -^(44)'iM. (1„ 

The independence of the expression (10) of the function A2 is a direct consequence of the gauge invariance. Then, in 
the one-loop approximation, 

W{C^; 6, MVA^ a,(M)) = 1 - ^^Cphix) ( J - 7b + ln47r + In ^) , (13) 

and the cusp dependent renormalization constant [23] , within the modified minimal subtraction scheme, reads 

Zcusp [C^; s, a,{n)] = 1 + ^^Cphix) - 7£ + ln47r) . (14) 

The detailed description of the renormalization procedure within the present approach has been made in [7] and will 
be omitted here for brevity. 

Using Eq. (10), one finds the known one-loop result for the perturbative field, which contains the dependence on 
the UV normalization point /i^ and IR cutoff A^ {e.g., [21]): 

<\Cx) = -^CW^(x)ln^. (15) 

Therefore, in the leading order the kinematic dependence of the expression (7) is factorized into the function h{x}, 
which at large-Q^ is approximated by 

/i(x)cxln^. (16) 

In this regime, the dependence of W on the UV normalization scale fj, (which can also be treated as an arbitrary 
factorization scale dividing the hard and soft subprocesses [24]) is governed by the renormalization group (RG) 
equation 



d ^. .d\d\n W{Q^) ^ .... 
% + ^^'^d-g ) d InQ^ = ["^^'^^J ' ^''^ 



where Tcuspics) is the universal cusp anomalous dimension evaluated in the perturbation theory. In Eq. (17), we 
take the logarithmic derivative in in order to avoid problems with light-cone singularities at = [24]. The 
solution of the RG equation leads to the evolution equation 



d lng2 
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where the function Wnp gives the initial condition at /i^ = and has to be found by the nonperturbative methods 
[25, 26]. Solving Eq. (18), we take the arbitrary upper bound for the squared momenta of soft gluons equal to the 
hard scale fj,^ = Q'^ and find 



In 



r'^' dx r r 

JQI X [j), 



d Wnp{x) 

d \nx 



(19) 



which immediately leads to the conclusion that the leading large-Q^ behavior of the quark form factor Fq{Q^) including 
all logarithmic corrections is controlled by the universal cusp anomalous dimension (17) and can be expressed in the 
following form (for comparison, see [24]): 



exp 



. r § ln^re„.p[a,(0]-lng § T^uspiasiO] +Wnp{Q 



(20) 



where Wq = W(Qf^) contains both pcrturbativc and nonperturbative contributions. From the one-loop result (15), 
the cusp anomalous dimension which satisfies the RG equation (17) in one-loop order is given by: 



r« (a.(M)) = . 



(21) 



Substituting the anomalous dimension (21) in the one-loop approximation for the strong coupling into the Eq. (20), 
one finds 



= exp 



^n^ln^''^^'/'^'^ 



00 A2 In(g2/A2) Ql \^ ln(AVA2) ' ' + 



(22) 



where A is the QCD scale. The singularity in Eq. (22) originates from the region where the IR cutoff approaches A, 
i. e., where the coupling constant increases, and then may have a nonperturbative nature. 

III. EFFECTS OF THE IR RENORMALONS 

In order to determine the structure of the nonperturbative function Wnp in Eqs. (20, 22), it is instructive to study 
the corrections due to IR renormalons [9] . In the present situation, one can expect the corrections proportional to the 
powers of both scales /z^ and A^. However, taking into account the evolution in /z^ to the hard characteristic scale of 
the process (19), we treat the power fi"^ terms to be strongly suppressed, and focus on the power A^ corrections. 
To find them, let us consider the perturbative function Ai(e, 0, /k^/A^) in the Eq. (10). The insertion of the fermion 
bubble 1-chain to the one-loop order expression (7) is equivalent to replacement of the frozen coupling constant g'^ by 
the running one g'^ — > g'^{k'^) = 47ras(fc^) [25]: 



Ai(£,0,/i7A^) = -47rK 



2e 



d'^k 
(27r)" 



e*'=^(5(z2) 
P(/e2 + A2) ■ 



(23) 



For the sake of convenience, we work here in Euclidean space. Using the integral representation for the one-loop 
running coupling asik"^) = d<T{h^ /k'^Y^ , b = I3q/Ait, we find 



Ai(£,0,/x7A') 



1 / ^£X f°° r(i-x-£)r(i + x + £) /a^ 



(24) 



To define properly the integral on the right hand sided of Eq. (24), one needs to specify a prescription to go around 
the poles, which are at the points S N. Of course, the result of integration will depend on this prescription 

giving an ambiguity proportional to (A^/A^) for each pole. Then, the IR renormalons produce the power corrections 
to the one-loop perturbative result, which we assume to exponentiate with the latter [25, 26]. Extracting from Eq. 
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(24) the UV singular part in vicinity of the origin a; = 0, wc divide the integration interval [0, oo] in two parts [0, S] 
and [(5, oo], where S < 1. This procedure allows us to evaluate the ultraviolet and the renormalon- induced pieces 
separately. For the ultraviolet piece, we apply the expansion of the integrand in Ai in powers of small x and replace 
the ratio of F functions by exp(— 7£;£): 



Ar(£,0,M7A^) 



/?o(l - e) 



E(-) 

k,n=0 



^(ln47r-7B + ln^ 



f 

Jo 



dx 



which after subtraction of the poles in the modified minimal subtraction scheme becomes 



Ar(0,M7A^) = 



In analogy with results of Ref. [27], this expression may be rewritten in a closed form as 



dx 



Substituting 



into Eq. (17) one finds 



/3o(l-e)yo ^ 
d VF(i)(Q2) 



-X In 



-X In - 



d lnQ2 

d lnW^(i)(Q') 
d lng2 



2Cir(l-e)Ar(0,/x7A^) 



= -FW (a,(^))(l-exp 
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(25) 



(26) 



(27) 



(28) 



(29) 



The second exponent in the last equation yields the power suppressed terms (A^/Q^) in large-Q^ regime. In the 
leading logarithmic approximation (LLA) Eq. (28) is reduced to 



dW^^^Q^) 
d lng2 



2Cf 



In 



ln(MVA2) 
ln(AVA2) 



(30) 



The last expression obviously satisfies the perturbative evolution equation (17). 

The remaining integral in Eq. (24) over the interval [S, oo] is evaluated at e = since there are no UV singularities. 
The resulting expression does not depend on the normalization point ii, and thus it is determined by the IR region 
including nonperturbative effects. It contains the renormalon ambiguities due to different prescriptions in going 
around the poles S„ in the Borel plane which yields the power corrections to the quark form factor. 

After the substitution fj,^ = and integration over d(ln Q^), we find in LLA (for comparison, see Eq. (22)) 



= exp 



2Cf 



ln(QVA2) 
A2 In(g2/A2) 



1 



1-ln 



ln(Q^/A2) 
ln(A2/A2) 



ln^0.e„(A2,A2) 



\Ql) 



(31) 



where the function (/)reri(A^, A2) — X)fc=o '^^^(^^/A^)'^ accumulates the effects of the IR renormalons, as well as the 
other nonperturbative information. The coefficients (j)k cannot be calculated in perturbation theory and can be 
treated as the minimal set of nonperturbative parameters. It is worth noting that the logarithmic dependence of 
the renormalon induced corrections in the large-Q2 regime is factorized, and thus the Eq. (31) reproduces exactly the 
structure of nonperturbative contributions found in the one-loop evolution equation (22) with respect to the large-Q^ 
asymptotic behavior. 

IV. LARGE-Q^ BEHAVIOUR OF THE INSTANTON INDUCED CONTRIBUTION 



Let us consider the instanton induced corrections to the perturbative result. The instanton field is given by 

fT° 1 6 

A^{x; p) = Al{x; p)^ = -R"''a"r;±^,(a; - zq)Mx - zq; p), 
^ 9 



(32) 
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where R"'' is the color orientation matrix [a = 1, (N^ — l),b = 1, 2, 3] , cr°'s are the Pauli matrices, and (±) corre- 
sponds to the instanton, or anti-instanton. The averaging of the Wilson operator over the nonpcrturbativc vacuum is 
reduced to the integration over the coordinate of the instanton center zq, the color orientation and the instanton size 
p. The measure for the averaging over the instanton ensemble reads dl = dR d'^Zo dn{p), where dR refers to the av- 
eraging over color orientation and dn{p) depends on the choice of the instanton size distribution. Taking into account 
Eq. (32), we write the Wilson integral (4), which defines the instanton induced contribution to the nonperturbative 
part in Eq. (20), in the single instanton approximation in the form 



where 



Wj(C^) = — (0|Trexp(ia»r)|0) , 



(33) 



(34) 



We omit the path ordering operator V in Eq. (33) because the instanton field (32) is a hedgehog in color space, and 
so it locks the color orientation by space coordinates. Although in certain situations, the integrals of this type (Eq. 
(34)) can be evaluated explicitly [3], the calculation of the total integral (34) for a given contour requires an additional 
work, so we must restrict ourselves with the weak- field approximation. In contrast to our previous paper [7], here we 
use the cutoff to regularize the IR divergences in the instanton case, while the UV divergences do not appear at 
all due to the finite instanton size. Then, in case of the instanton field, the LO contribution reads 



W^'\C^) = 2h{x) J dn{p) Ai{0,p'X') , 



where 



A((0,p^A^) 



d'^k 
(2^ 



e''''S{z^) [2<^'(fc2;p) 



(35) 



(36) 



Here, (p{k'^\p) is the Fourier transform of the instanton profile function ip{z'^;p) and (f'{k^;p) is it's derivative with 
respect to fc^. Note, that for the instanton calculations, it is necessary to map the scattering angle x to the Euclidean 
space by the analytical continuation [28] x — > 17 , and perform the inverse transition to the Minkowski space-time in 
the final expressions in order to restore the dependence. In the singular gauge, when the profile function is 



one gets 



where 



and 



A( (0, p'X') = ^ [Hp'X') $o(p^A2) + Mp'X')] , 



dz 



2 \ 2 2 

1 + e'' ^ - 2e^-'' 



$i(p2x2) = V / dxdydz 



z{l - z) 

[-p^X^jxz + y{l - ^))]" p^A^[^^+„(i-^)] 



lim $o(p A^) = 1 



n\n 



lim $i(p2A2) =0 



are the IR-finite expressions. At high energy the instanton induced contribution is reduced to the form 

dwi{Q^) TT^ 



^ J dn{p) [Hp'X') $o(p'A2) + $i(p'A2)] = -Bi{X') . 



d lnQ2 

Here we used the exponentiation of the single- instanton result in a dilute instanton ensemble [7]: 



(37) 



(38) 



(39) 



(40) 



(41) 



Wi = exp (wi) , 



(42) 
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and took only the LO term of the weak-field expansion (7): W^^'' wj. 

In order to estimate the magnitude of the instanton induced effect we consider the standard distribution function 
[29] supplied with the exponential suppressing factor, which has been suggested in Rcf. [30] (and discussed in Ref. 
[31] in the framework of constrained instanton model) in order to describe the lattice data [32] 



dn{p) = ^ Cn, 



2tt 



2Na 



exp 



27r 



where the constant is 



Cm 



0.466 e-i-679^- 
(A^e-l)!(iVe-2)! 



0.0015 , 



(43) 



(44) 



(7 is the string tension, f3 = Pq + O [Q!s(/i, )], and /i,. is the normalization point [34]. Given the distribution (43) the 
main parameters of the instanton liquid model — the mean instanton size p and the instanton density n — will read 



r(/3/2 - 3/2) 1 
r(/3/2-2) V2^ 



(45) 



n = 



Cn,T{P/2 -2) ( 27r 



Q!s(/9-l) 



2JVc 



A, 



12^ ) 



(46) 



In Eq. (46) wc choose, for convenience, the normalization scale /i,- of order of the instanton inverse mean size p^^ . 
Note, that these quantities correspond to the mean size pq and density no of instantons used in the model [33] , where 
the size distribution (43) is approximated by the delta function dn(p) = no6{p — po)dp . 
Thus, we find the leading instanton contribution (41) in the form 



Bi = Kn^fip^ In 



27ro- 
1^ 



1 



\2T:a J 



where 



K = 



r(/3o/2)[r(/3o/2-2)]= 



2 [r(/3o/2 - 3/2)]4 
and we used the one loop expression for the running coupling constant 

2n „ UN, 



0.74 



Oisip = - 



/3oln pA 



/3o 



2nt 



(47) 



(48) 



(49) 



The packing fraction n^np'^ characterizes diluteness of the instanton liquid and within the conventional picture its 
value is estimated to be 0.12 , if one takes the model parameters as (see Ref. [10]) 



1 fm ^, p^ 1/3 fm 



(0.44 GeVf 



(50) 



The leading contribution to the quark form factor at asymptotically large is provided by the (perturbative) 
evolution governed by the cusp anomalous dimension (21). Thus, the instantons yield subleading effects to the 
large-Q^ behavior accompanied by a numerically small factor 



Bi w 0.02 , 



(51) 



as compared to the perturbative term 2Cf/(3o « 0.24. 

Therefore, from Eqs. (41) and (31), we find the expression for the quark form factor at large Q"^ with the one-loop 
perturbative contribution and the nonperturbative contributions (the function W„p in Eq. (22)) which include both 
the IR renormalon and the instanton induced terms 



exp 



2Cf 



ln(QVA2) 



1 



1 - In 



ln(A2/A2) 



91 



(52) 
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It is clear, that while the asymptotic (double-logarithmic) behavior is controlled by the perturbative cusp anomalous 
dimension, the leading nonperturbative corrections results in a finite renormalization of the next-to-leading (logarith- 
mic) perturbative term. From the formal point of view, the evolution equation (20) describing the large-Q^ asymptotic 
is valid even at the low scales ~ 1 GeV^, since the only condition of applicability of the Wilson integrals approach 
is >> }?,m?. However, in the low-energy domain the perturbative one- loop cusp anomalous dimension Tcusp{ois) 
(21) should be supplemented by higher loop corrections, and thus the explicit formula (52) would include additional 
logarithmic terms. The relevance of the instanton induced part (41) in the low-energy domain calculated in the dilute 
gas approximation can be questioned, and the additional consideration within more proper framework may help one 
to verify it. Indeed, the corrections to the single instanton approximation may be large at sufficiently low momenta 
(see for recent discussions, e.g. Ref. [35] and references therein). At the moment, we can say confidently that the 
evolution equation is valid at > 1 GeV'^. 

We have to comment that the weak field limit used in the instanton calculations may deviate from the exact 
result. Nevertheless, we expect that using of the instanton solution in the singular gauge, that concentrate the field 
at small distances, leads to the reasonable numerical estimate of the full effect. Thus, the resulting diminishing of 
the instanton contributions with respect to the perturbative result appears to be reasonable output. It should be 
emphasized that in the present paper, all the calculations have been performed analytically while the evaluation 
of the instanton contribution beyond the weak field approximation requires a numerical analysis, which will be the 
subject of a separate work. Moreover, the use of the singular gauge for the instanton solution allows us to prove the 
exponentiation theorem for the Wilson loop in the instanton field [7] which permits one to express the full instanton 
contribution as the exponent of the all-order single instanton result (42). 

V. CONCLUSION 

We analyzed the structure of the nonperturbative corrections to the quark form factor at large momentum transfer. 

In order to model the nonperturbative effects, we studied the quark scattering process in the backgroimd of the 
instanton vacuum. The instanton induced contribution to the color singlet quark form factor is calculated in the 
large momentum transfer regime. It was shown that the instanton induced corrections correspond to the leading term 
proportional to In . The magnitude of these corrections is determined by the small instanton liquid packing fraction 
parameter, and they can be treated as finite renormalization of the subleading perturbative part (52). In addition 
to this, the minimal set of the nonperturbative parameters is found considering infrared renormalon ambiguities of 
the perturbative series. Within this approach, it is shown that the leading large-Q^ behavior of nonperturbative 
contributions should also be determined by the logarithmic term ~ In Q^, what is consistent with the instanton 
analysis. 

Let us emphasize that our results are quite sensitive to the prescription how to make the integration over instanton 
sizes finite. For example, if one used the sharp cutoff then the instanton would produce strongly suppressed power 
corrections such as cx {A/Q)^°. However, we think that the distribution function (43) should be considered as more 
realistic, since it reflects more properly the structure of the instanton ensemble modeling the QCD vacuum. Indeed, 
this shape of distribution was recently advocated in Refs. [30, 31] and supported by the lattice calculations [32]. 

Finally, we think that the instanton induced effects are more interesting for theoretical investigation and more 
important for phenomenology in the hadronic processes which possess two energy scales, such as the total center-of- 
mass energy s (hard characteristic scale), and the squared momentum transfer —t which is small compared to the 
latter: —t« s. One of the most interesting examples of such processes is the parton-parton scattering and the soft 
Pomeron problem [3, 4, 11]. Another important situations where the nonperturbative (including instanton induced) 
effects can emerge are the transverse momentum distribution of vector bosons in the Drell-Yan process (see, e.g., Refs. 
[25, 26]), and the phenomenon of saturation in deep-inelastic scattering (DIS) at small x [36]. The explicit evaluation 
of the instanton effects in these processes will be the subject of our forthcoming study. 
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